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Modular structure is ubiquitous among real-world networks from related proteins to social groups. 
Here we analyze the modular organization of brain networks at a large-scale (voxel level) extracted 
from functional magnetic resonance imaging (fMRI) signals. By using a random walk-based method, 
we unveil the modularity of brain-webs, and show modules with a spatial distribution that matches 
anatomical structures with functional significance. The functional role of each node in the network 
is studied by analyzing its patterns of inter- and intra-modular connections. Results suggest that the 
modular architecture constitutes the structural basis for the coexistence of functional integration of 
distant and specialized brain areas during normal brain activities at rest. 



PACS numbers: 89.75.-k, 87.19.lf, 87.19.1m 

There is a growing interest in studying the con- 
nectivity patterns extracted from brain signals 
during different mental states. Current stud- 
ies suggest that brain architecture leads neural 
assemblies to be coordinated with an optimized 
wiring cost. Brain webs coordinate a mosaic of 
brain modules, carrying out specific functional 
tasks and integrated into a coherent process. We 
analyze the modular structure of brain networks 
extracted from fMRI signals in humans at rest. 
Using a random walk-based method we identify 
a non-random modular architecture of brain con- 
nectivity. This approach is fully data driven and 
relies on no a priori choice of a seed brain region 
or signal averaging in predefined brain areas. The 
analysis of intra- and inter-modules connections 
leads us to relate a node's connectivity to a lo- 
cal information processing, or to the integration 
of distant anatomo/functional brain regions. We 
also find that the spatial distribution of the re- 
trieved modules matches with brain areas associ- 
ated with specific functions, assessing a functional 
significance to the modules. In our conclusions, 
we argue that a modular characterization of the 
functional brain webs constitutes an interesting 
model for the study of brain connectivity during 
different pathological or cognitive states. 



I. INTRODUCTION 

From the brain over the Internet to social groups, com- 
plex networks are a prominent framework to describe col- 
lective behaviors in many areas p] . Many of real- world 
networks exhibit topological features that can be cap- 
tured neither by regular connectivity models as lattices, 
nor by random configurations \2j 2^ ■ Under this frame- 
work, recent studies of complex brain networks have at- 
tempted to characterize the connectivity patterns ob- 



served under functional brain states. Electroencephalog- 
raphy (EEG), magnetoencephalography (MEG), or func- 
tional magnetic resonance imaging (fMRI) studies have 
consistently shown that human brain functional networks 
during different pathological and cognitive neurodynam- 
ical states display small world (SW) attributes [HO [SI 13 
IHJ [5] • SW networks are characterized by a small average 
distance between any two nodes while keeping a rela- 
tively highly clustered structure. Thus, SW architecture 
is an attractive model for brain connectivity because it 
leads distributed neural assemblies to be integrated into 
a coherent process with an optimized wiring cost [10| [11] . 

Another property observed in many networks is the 
existence of a modular organization in the wiring struc- 
ture. Examples range from RNA structures to biological 
organisms and social groups. A module is currently de- 
fined as a subset of units within a network such that con- 
nections between them are denser than connections with 
the rest of the network. It is generally acknowledged that 
modularity increases robustness, flexibility and stability 
of biological systems [HI [T5] . The widespread character 
of modular architecture in real-world networks suggests 
that a network's function is strongly ruled by the orga- 
nization of their structural subgroups. 

Empirical studies have lead to the hypothesis that spe- 
cialized neural populations are largely distributed and 
linked to form a web-like structure [13]. The emergence 
of any unified brain process relies on the coordination of 
a scattered mosaic of modules, representing functional 
units, separable from -but related to- other modules. 
Characterizing the modular structure of the brain may 
be crucial to understand its organization during different 
pathological or cognitive states. 

Previous studies over the mammalian and human brain 
networks have successfully used different methods to 
identify clusters of brain activities. Some classical ap- 
proaches, such as those based on principal components 
analysis (PCA) and independent components analysis 
(ICA), make very strong statistical assumptions (or- 
thogonality and statistical independence of the retrieved 
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components, respectively) with no physiological justifi- 
cation [m [16]. Although a number of studies investi- 
gating the organization of anatomic and functional brain 
networks have shown very interesting properties of the 
macro-scale brain architecture [T71 \TE\, little is known 
about the network structure at a finer scale (at a voxel 
level) . Current approaches are based on the use of a pri- 
ori coarse parcellations of the cortex [IJ IH] ; or on partial 
networks defined by a seed voxel [H]. Nevertheless, seed- 
based descriptions may fail to describe the global behav- 
ior of the brain, as they only consider the connectivity 
of the reference voxel. On the other hand, parcellation 
schemes reduce the analysis to a macro-scale fixed by an 
a priori definition of the brain areas. Further, a recent 
study shows that the topological organization of brain 
networks is affected by the different parcellation strate- 
gies applied [20] . 

Here we focus on a completely data-driven framework 
to study the connectivity of brain networks extracted 
directly from functional magnetic resonance imaging 
(fMRI) signals at voxel resolution. A random walk-based 
algorithm is used to assess the modular organization of 
functional networks from healthy subjects in a resting- 
state condition. Results reveal that functional brain webs 
present a large-scale modular organization significatively 
different from that arising from random configurations. 
Further, the spatial distribution of some modules fits well 
with previously defined anatomo-functional brain areas, 
assessing a functional significance to the retrieved mod- 
ules. Based on the patterns of inter- and intra-modular 
connectivities, we also study the roles played by differ- 
ent brain sites [211 . Results provide a characterization of 
the functional scaffold that underly the coordination of 
specialized brain systems during spontaneous brain be- 
havior. 



II. DATA ADQUISITION AND 
PREPROCESSING 



BOLD fMRI data were acquired using a T2*-weighted 
imaging sequence during a period of 10 minutes from 7 
healthy right-handed subjects. The study was performed 
with written consent of the subjects and with the ap- 
proval of local ethics committees. During the scan, all 
subjects were instructed to rest quietly, but alert, and 
keep their eyes closed. 500 volumes of gradient echopla- 
nar imaging (EPI) data depicting BOLD contrast were 
acquired. In the acquisition, we used the following pa- 
rameters: number of slices, 21 (interleaved); slice thick- 
ness, 4 mm; inter-slice gap, 1 mm; matrix size, 64 x 64; 
flip angle, 90 °; repetition time (TR), 1250 ms; echo time, 
30 ms; in-plane resolution, 3x3 mm^. Subsequently, a 
high resolution structural volume was acquired via a Tl- 
weightcd sequence (axial; matrix 192 x 256 x 160; FOV 
192 X 256 X 160 mm^; slice thickness; 1 mm; in-plane 
voxel size, 1 x 1 mm^; flip angle 15° ; TR, 1620 ms, TI, 



950 ms; TE, 3.87 ms) to provide the anatomical reference 
for the functional scan. 

All acquired brain volumes were corrected for mo- 
tion and differences in slice acquisition times using the 
SPM5 (http://www.fil.ion.ucl.ac.uk) software package. 
After correction, fMRI datasets were coregistered to the 
anatomical dataset and normalized to the standard tem- 
plate MNI, enabling comparisons between subjects. Due 
to computational limitations, normalized and corrected 
functional scans were subsampled to a 4x4x4 mm res- 
olution, yielding a total of 20898 voxels (nodes in the 
network). To eliminate low frequency noise (e.g. slow 
scanner drifts) and higher frequency artifacts from car- 
diac and respiratory oscillations, time-series were digi- 
tally filtered with a finite impulse response (FIR) filter 
with zero-phase distortion (bandwidth 0.01 — 0.1 Hz) [19]. 



III. ESTIMATION OF FUNCTIONAL 
CONNECTIVITY 



A functional link between two time series Xi(t) and 
Xj{t) (normalized to zero mean and unit variance) was 
defined by means of the linear cross-correlation coeffi- 
cient computed 8jS Tj^j - — - (^X^ {t)xj{t)), where (•) denotes 
the temporal average. For the sake of simplicity, we only 
considered here correlations at lag zero. To determine 
the probability that correlation values are significantly 
higher than what is expected from independent time se- 
ries, rij(O) values (denoted r^) were firstly transformed 
by the Fisher's Z transform 



Z„ =0.51n( 

1 - r,. 



(1) 



Under the hypothesis of independence, Zij has a nor- 
mal distribution with expected value and variance 
l/{df — 3), where df is the effective number of degrees 
of freedom [211 HSl HI] . If time series are formed of in- 
dependent measurements, df simply equals the sample 
size, N. Nevertheless, autocorrelated time series do not 
meet the assumption of independence required by the 
standard significance test, yielding a greater Type I er- 
ror [55] US] |M] . In presence of auto-correlated time series 
df must be corrected by the following approximation: 
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where r^xiT) is the autocorrelation of signal x at lag 
r. Other estimators of df, and statistical significance 
tests for auto-correlated time series can be found in [25] . 
To correct for multiple testing, the False Discovery Rate 
(FDR) method was applied to each matrix of r^j val- 
ues [26 . With this approach, the threshold of signifi- 
cance rth was set such that the expected fraction of false 
positives is restricted to g < 0.001. 

In the construction of the networks, a functional con- 
nection between two brain sites was assumed as an undi- 
rected and unweighted edge (Aij = 1 if rij > rth; 
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FIG. 1: Cumulative degree distributions P{K > k) estimated 
from all subjects. The inset depicts the average degree distri- 
bution. Black and gray solid lines correspond to the observed 
distribution and the fitted truncated power law model, re- 
spectively. Dotted, dashed and dot-dash lines correspond to 
fitted power law, exponential law and truncated Pareto law 
distributions, respectively. 

and zero otherwise). Although topological features can 
also be straightforwardly generalized to weighted net- 
works, we obtained qualitative similar results (not re- 
ported here) for weighted networks with a functional con- 
nectivity strength between nodes given by Wij = Vij . 

To characterize the topological properties of a network, 
a number of parameters have been described. Here we 
use three key parameters: mean degree {K), clustering 
index C and global efficiency £' [U El |3] . Briefly, the de- 
gree ki of node i denotes the number of functional links 
incident with the node and the mean degree is obtained 
by averaging ki across all nodes of the network. The clus- 
tering index quantifies the local density of connections in 
a node's neighborhood. For a node i, the clustering co- 
efficient Ci is calculated as the number of links between 
the node's neighbors divided by all of their possible con- 
nections and C is defined as the average of Ci taken over 
all nodes of the network. The global efficiency E pro- 
vides a measure of the network's capability for informa- 
tion transfer between nodes and is defined as the inverse 
of the harmonic mean of the shortest path length L^j 
between each pair of nodes. 

Figure [T] shows superimposed the degree distribu- 
tions for the seven studied subjects. For each net- 
work, goodness-of-fit was compared here using Maxi- 
mum Likelihood methods and the Kolmogorov-Smirnov 
statistic (KS) for four possible forms of degree distri- 
bution p{k): a power law p{k) oc k an exponential 
p(k) oc exp"'*''^; a truncated Pareto p{k) oc — 
(^a+i^-i^a. g^jjj g^jj exponentially truncated power law 
p{k) oc exp(— fc/fcc). The bestfitting were obtained 
for the truncated power law (KS = 0.0421 compared with 
KS = 0.1028, 0.2632 and 0.3278 for the exponential law, 
the truncated Pareto and the power law distribution, 
respectively). Estimated parameters for the truncated 
power law are a = 0.7688 ± 0.1455, kc = 410 ± 351. 

Values of the topological parameters are summarized 
in Table |l] To asses the statistical significance of brain 
connectivity, we perform a benchmark comparison of the 
functional connectivity patterns. For this, the topologi- 
cal features of brain webs are compared with those ob- 
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710.65 248.37 815.17 263.59 134.69 133.06 201.16 
0.4954 0.3901 0.4865 0.3856 0.3541 0.3389 0.3638 
0.0340 0.0119 0.0390 0.0126 0.0064 0.0064 0.0096 
0.3888 0.3569 0.4135 0.3447 0.3104 0.3132 0.3269 
0.5170 0.4973 0.5195 0.5004 0.4337 0.4322 0.4810 



TABLE I: Parameters for real and randomized networks: {K), 
mean degree; C, clustering index; E, global efficiency; Omd 
denotes the average of parameter 6 obtained from 10 equiva- 
lent randomized networks. Single asterisks indicate that this 
parameter has a significance level of p < 10~*. 



taincd from equivalent random wirings. To create an 
ensemble of equivalent random networks we use the algo- 
rithm described in |T]. According to this procedure, each 
edge of the original network is randomly rewired avoiding 
self- and duplicate connections. The obtained random- 
ized networks thus preserve the same mean degree as the 
original network, whereas the rest of the wiring structure 
is random. The significance of a given topological param- 
eter 6 was assessed by quantifying its statistical deviation 
from values obtained for the ensemble of randomized net- 
works. Let /i and a be the mean and standard deviation 
of the parameter 9 computed from such an ensemble. The 
significance is given by the ratio S = (9 — 11) /a whose 
p-value is given by the Chebyshev's inequality (for any 
statistical distribution of 9: p{\T,\ ^ () ^ 1/C^ where C is 
the chosen statistical threshold) [27] . 

The sparse connectivity of functional brain networks 
was found to be significatively different from randomized 
wirings for all the subjects. Brain networks yielded larger 
clustering values {p < 10~^) than the equivalent random 
configurations, but similar efficiency values, indicating a 
connectivity with SW attributes. 



IV. MODULAR ANALYSIS OF BRAIN 
NETWORKS 

A potential modularity of brain-webs is suggested by 
the fact that brain networks display a clustering index 
approximately one order of magnitude larger than that 
obtained from random configurations [28 . Although the 
notion of module results very intuitive, in general it is 
difficult to define formally. It is currently accepted that 
a partition V — {Ci, . . . ,Cm} represents a good modu- 
lar structure if the portion of edges inside each module 
Ci (intra-modular edges) is high compared to the por- 
tion of edges between them (inter-modular edges) . The 
modularity Q{V), for a given partition V oi s. network is 
formally defined as [29] : 
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where M is the number of modules, L is the total num- 
ber of connections in the network, Ig is the number of 
connections between vertices in module s, and kg is the 
sum of the degrees of the vertices in module s. 
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Qrnd 


0.4385 
0.0065 


0.5814 0.4223 0.5538 0.5648 0.5749 
0.0169 0.0057 0.0160 0.0274 0.0279 


0.5362 
0.0201 



TABLE II: Modularity for real (Q) and randomized networks. 
Qrnd denotes the average obtained from 10 equivalent ran- 
domized networks. Double asterisks denotes a significance 
level of p < 10"^ 

To partition the functional networks in modules, we 
used a random walk-based algorithm [3D], because of its 
ability to manage very large networks, and its good per- 
formances in benchmark tests 30, 31,. In a nutshell, a 
random walker on a connected graph tends to remain into 
densely connected subsets corresponding to modules. Let 
Pij = to be the transition probability from node i to 
node j, where Aij denotes the adjacency matrix and ki is 
the degree of the i*^ node. This defines the transition ma- 
trix {P^)ij for a random walk process of length t (denoted 
here P/^ for simplicity) . The metric used to quantify the 
structural similarity between vertices is given by 
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Using matrix identities, the distance p can be written as 
Pij = YJa=2^'^ai'^a{i) - Va{j))^; where (AQ)i^a^„ and 
{va)i^a^n ^rc the n eigenvalues and right eigenvectors 
of the matrix P, respectively [3U] . This relates the ran- 
dom walk algorithm to current methods using spectral 
properties of the graphs [32^ I33|. The current approach, 
however, needs not to explicitly compute the eigenvectors 
of the matrix; a computation that rapidly becomes in- 
tractable when the size of the graphs exceeds some thou- 
sands of vertices. 

To find the modular structure, the algorithm starts 
with a partition in which each node in the network is the 
sole member of a module. Modules are then merged by 
an agglomerative approach based on a hierarchical clus- 
tering method. Following Ref. [30], if two modules Ci and 
C2 are merged into a new one C3 = Ci UC2, the transition 

matrix is updated as follows: F^^j, = — 

where \Ci\ denotes the number of elements in module Ci. 
The algorithm stops when all the nodes are grouped into 
a single component. At each step the algorithm evaluates 
the quality of partition Q. The partition that maximizes 
Q is considered as the partition that better captures the 
modular structure of the network. In the calculation of 
Q, the algorithm excludes small isolated groups of con- 
nected vertices without any links to the main network. 
However, these isolated modules are considered here as 



part of the network for the calculation of the topological 
parameters. 

As reported in Table [TTj a modular structure is con- 
firmed by the high values of Q obtained for the optimal 
partition of the networks (a value of Q > 0.3 is in prac- 
tice a good indicator of modularity in a network [M]V 
Further, values of modularity for all the subjects were 
statistically significant when compared with randomized 
wirings {p < 10^^). To assess the stability of the par- 
tition structure across subjects we used the Rand index 
J [35 , which is a traditional criterion for comparison of 
different results provided by classifiers and clustering al- 
gorithms, including partitions with different numbers of 
classes or clusters. For two partitions P and P' the Rand 
index is defined as J = — ttt^t^; where a is number of 

a-\-b-\-c+a ' 

pairs of data objects belonging to the same class in P 
and to the same class in P' , b is number of pairs of data 
objects belonging to the same class in P and to different 
classes in P', c is the number of pairs of data objects be- 
longing to different classes in P and to the same class in 
P', and d is number of pairs of data objects belonging to 
different classes in P and to different classes in P'. Thus 
index J yields a normalized value between (if the two 
partitions are randomly drawn) and 1 (for identical par- 
tition structures). For our data, the values of J indicate 
a moderate stability of the partition structure across all 
subjects (J = 0.5148). 

To assess a functionality to the different groups of the 
modular brain webs, we compared the spatial distribu- 
tion of the recovered modules with a previously reported 
anatomical parcellation of the human brain |36j . For 
the sake of simplicity, we only consider here communities 
whose size was larger than 40 voxels (^ 0.2% of the size 
of the whole network), which yields Nc — 22 modules. 

Fig. [2]illustrates the spatial distribution of the modules 
retrieved from the averaged connectivity matrix com- 
puted over all subjects. Results show that the spatial 
distribution of recovered modules fits well some brain 
systems. Module 22 for instance, includes 75% of the 
primary visual areas VI, while module 5 overlaps half of 
the ventral visual stream (brain areas V2 and V4), and 
visual areas of the V3 region (cuneus and precuneus) are 
included (~ 40%) in the module 4. Module 20 includes 
most of the subcortical structures caudate and thalamus 
nuclei (covered at 70% and 75%, respectively). The au- 
ditory system is included by module 12 that overlaps 
primary and secondary areas plus associative auditory 
cortex (60 — 70%). Modules 11, 16 and 21 cover most 
(40 — 70%) of the somatosensory and motor cortices; and 
language related areas are mainly included (> 60%) in 
module 10. 

Importantly, some modules include distant brain lo- 
cations that are functionally related, e.g. the language 
related areas (modules 10), the auditory system (module 
12), or brain regions involved in high level visual process- 
ing tasks (module 5). This spatially distributed organiza- 
tion of modules rules out the possibility that modularity 
simply emerges as a consequence of vascular processes 
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FIG. 2: Main functional brain modules: brain sites belonging to each module were coloured and superimposed onto an 
anatomical image. The sagittal anatomical images at the top right of each plot indicate the relative position of imaged slices 
of each row. For the sake of clarity, we show only those communities with a size larger than 40 voxels (~ 0.2% of the size of 
the whole network). 



or local physiological activities independent of neuronal 
functions [571f38j. 

Modules assignment provides the basis for the classifi- 
cation of nodes according to their patterns of intra- and 
inter-modules connections, which conveys significant in- 
formation about the importance of each node within the 
network [2T] . 

The within-module degree z-score measures how well 
connected the node i is to other nodes in the module, 
and is defined as: 



(5) 



where ki is the number of links of node i to other nodes 
in its module Si, kg. is the average of k over all the nodes 
in Si , and ak^ is the standard deviation of A; in . Thus 
node i will display a large value of Zi if it has a large 
number of intra-modular connections relative to other 
nodes in the same module, i.e. it measures how well 
connected a node is to other nodes in the module). 



The extent a node i connects to different modules is 
measured by the participation coefficient pci defined as: 



PCi 



M 

E 



(6) 



where fcj^ is the number of links of node i to nodes in 
module s, and ki is the degree of node i. The participa- 
tion coefficient takes values of zero if a node has most of 
its connections exclusively with other nodes of its mod- 
ule. In contrast, pCi ^ 1 if their links are distributed 
among different modules in the network. 

The role (R^) of a node in the network can be assessed 
by its within-module degree and its participation coef- 
ficient, which define how the node is positioned in its 
own module and with respect to other modules 21J. Fig- 
ure [3] shows the distribution of the roles obtained from 
all the analyzed networks over the z — pc parameter 
space. Most of the nodes in the functional brain networks 
(~ 98%) can be classified as non-hubs (indicated by the 
gray area in Fig. [3j(b)), while only a minority of them 
are module hubs 2%). Non-hubs nodes were classified 
as ultra-peripheral (Rl, 10.33%) having all their links 
within their own modules; peripherals (R2, 73.49%) with 
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FIG. 3: Role determination, as represented in the z — pc pa- 
rameter space, (a) Average density landscape (computed over 
all networks) depicted in logarithmic scale, (b) Histograms 
and error bars corresponding to the proportion of nodes for 
each role. Histograms are coloured according to the roles de- 
picted in the z — pc parameter space. 

most links within their modules; or non hub-connectors 
(R3, 13.67%) with half of their links to other modules. 
This distribution of roles strongly contrasts with that 
obtained from random configurations (results not show) 
where most nodes have their links homogeneously dis- 
tributed among all modules (R4 and R7). 

The anatomical distribution of the parameters z and 
pc is depicted in Figure |4j Interestingly, this represen- 
tation shows that the wiring structure of the brain has 
a non-homogeneous organization in terms of the z — pc 
parameters distribution. Examples of the different be- 
haviours that can be observed are: i) subcortical struc- 
tures (indicated by the orange arrow) display relatively 
high values for both z-score and pc parameters, indicat- 
ing a dense inter- and intra- modular connectivity ; ii) 
nodes belonging to brain areas associated to the primary 
visual system (pointed by the red arrow) have a scat- 
ter connectivity, yielding low values for both pc and z 
parameters; in) precuneus and cyngular gyrus areas (in- 
dicated by the yellow arrow) have a dense intra-modular 
connectivity (high values of z ) but few links to other 
modules (low values of pc); iv) frontal areas and some vi- 
sual regions related to associative functions (cian arrow) 
present more connections to other modules, which is re- 
flected in their low values of z and relatively high values 
of pc. 

V. CONCLUSION 

In conclusion, here we address a fundamental prob- 
lem in brain networks research: whether the sponta- 
neous brain behavior relies on the coordination (integra- 
tion) of a complex mosaic of functional brain modules 
(segregation). By using a random walk-based method 
we have identified a non-random modular structure of 
functional brain networks. In contrast to current ap- 
proaches m El , our procedure requires neither of sig- 
nal averaging in predefined brain areas, nor the defini- 
tion of seed regions, nor subjective thresholds to assess 
the connectivities. To our knowledge, this work provides 
the first evidence of a modular architecture in functional 
human brain networks at a voxel level. 



The modularity analysis of large-scale brain networks 
unveiled a modular structure in the functional connec- 
tivity. Although a one-to-one assignment of anatomi- 
cal brain regions to each detected module is difficult to 
define, results reveal a strong correlation between the 
spatial distribution of the modules and some well-known 
functional systems of the brain, including some of the fre- 
quently reported circuits underlying the functional activ- 
ity at rest [19 . It is worth to notice that, although the 
functional brain connectivity is strongly shaped by the 
underlying anatomical wiring (e.g. by the white matter 
pathways), future studied are needed to clearly examine 
the interplay between the structural substrate and the 
modular connectivity inferred from brain dynamics |39j . 

Our findings are in full agreement with previous studies 
about the structure of human brain networks. First, we 
have confirmed the degree distribution presents a power- 
law behavior over a wide range of scales, implying that 
there are a small number of regions with a large number 
of connections. We also found that brain connectivity 
shows a degree of clustering that is one order of mag- 
nitude higher than that of the equivalent random net- 
works while keeping similar efficiency values, suggesting 
that spontaneous brain behavior involves an optimized 
(in a SW sense) functional integration of distant brain 
regions [IIISJIS]. Further, the intrinsic non-random mod- 
ular structure suggested by the high values of the clus- 
tering index of brain networks was confirmed by a high 
degree of modularity obtained for the ensemble of sub- 
jects. 

Although the mechanisms by which modularity 
emerges in complex networks are not well understood, 
it is widely believed that the modular structure of com- 
plex networks plays a critical role in their functional- 
ity [211 30] • Functional brain modules can be related 
to a local -segregate- information processing while inter- 
modular connections allows the integration of distant 
anatomo/functional brain regions [41' . On the other 
hand, the SW and scale-free characteristics of brain webs 
provide an optimal organization for the stability, ro- 
bustnes, and transfer of information in the brain [6l[7l[8]. 
The modular structure constitutes therefore an attractive 
model for the brain organization as it supports the coex- 
istence of a functional segregation of distant specialized 
areas and their integration during spontaneous brain ac- 
tivity [351133]. Although the study of anatomical brain 
networks is a current subject of research, we suggest that 
a modular description might provide new insights into 
the understanding of human brain connectivity during 
pathological or cognitive states. 
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FIG. 4: Anatomical distribution of (left) the averaged within-module degree z, and (right) the averaged participation coefficient 
pc indices computed across all subjects. The sagittal anatomical images at the top right of each plot indicate the relative position 
of imaged slices of each row. 
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